Semifluxon molecule under control 
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Josephson junctions with a phase drop n in the ground state allow to create vortices of supercurrent carrying 
only half of the magnetic flux quantum $o ~ 2.07 x 10~^^ Wb. Such semifluxons have two-fold degenerate 
ground states denoted | (with flux +$0/2 and supercurrent circulating clockwise) and j (with flux — $o/2 and 
supercurrent circulating counterclockwise). We investigate a molecule consisting of two coupled semifluxons in 
a O-TT-O long Josephson junction. The fluxes (polarities) of semifluxons are measured by two on-chip SQUIDs. 
By varying the dc bias current applied to the O-tt-O junction, we demonstrate controllable manipulation and 
switching between two states, || and ||, of a semifluxon molecule. These results provide a major step towards 
employing semifluxons as bits or qubits for classical and quantum digital electronics. 

PACS numbers: 74.50.+r, 85.25.Cp 05.45.-a 74.20.Rp 



Introduction. Any introductory book on superconductiv- 
ity claims that magnetic flux in a superconducting loop is 
quantized with the single flux quantum being $0 ~ 2.07 x 
10~i^ Wb. The flux quantization owes to the 2tt periodicity 
of the phase of the wave function of the superconducting con- 
densate and leads to the existence of fluxons — Abrikosov 
or Josephson vortices of supercurrent carrying ±^q. In this 
case the phase grows smoothly by ±27r if one goes around a 
loop containing the vortex. It turns out that using a so-called 
TT Josephson junction (JJ), i.e. a JJ having a phase drop of tt 
in the ground state, embedded in this loop, a phase drop tt can 
be "stolen" (squeezed inside the tt junction) and the current 
circulating in a loop should provide only another tt or — tt to 
have a phase change equal to a multiple of 27r when one goes 
around the vortex. As a result, the magnetic flux threading the 
loop is equal to ±$0/2 — a semifluxon or antisemifluxon, 
respectively. The investigation of semifluxons intensified af- 
ter the advent of O-tt JJs that can be fabricated using d-wave 
superconductors m 13, [stlj using conventional superconductors 
with ferromagnetic barriersjJIst] or even creating an artifi- 
cial discontinuity of the Josephson phase using a pair of tiny 
current injectors attached to one of the electrodes of a conven- 
tional J J with insulating baiTier@01. 

A O-TT JJ consists of two connected parts: a part and a 
TT part. The part, if taken separately, would have a phase 
drop i^L ~ in the ground state, while the tt part would have 
/I = TT. If and TT parts are connected, the ground state phase 
should "decide" which value to take. If both and tt parts are 
long enough in terms of the Josephson penetration depth Aj 
(typically A j ^ 1 . . . 200 /im), the ground state phase ii{x) 
smoothly changes from to zLtt in the Aj-vicinity of a 0- 
TT boundary 013, fioll . This results in the appearance of a 
magnetic field cx d^/dx localized in the vicinity of the O-tt 
boundary and a supercurrent ^ sin ii{x) circulating around 
the O-TT boundary. The total magnetic flux associated with 
the localized field (created by a Josephson vortex of supercur- 
rent) is equal to $ = ±<i>o/2 (if the J J length L > Aj). The 
± sign stands for a doubly degenerate ground state with the 
supercurrent circulating clockwise or counterclockwise ifioll . 



Since semifluxons are pinned at the O-tt boundary and may 
have two opposite polarities, which we denote as | and j, 
they can be used in digital electronics to store and pro- 
cess information as classical or quantum bits. The physics 
of single semifluxons has been addressed by several groups 
theoretically S H EJll El El [3 H [3 [B H El 
and experimentally ||4|,|5l|7|]. However, in the systems of sev- 
eral coupled semifluxons, i.e. semifluxon molecules, only the 
ground states were imaged iB Eoll . Particularly interesting is 
the simplest system of two coupled semifluxons in a O-tt-O JJ, 
which has two degenerate ground states f j and jf. This ar- 
rangement was proposed as a candidate for semifluxon based 



bits and qubits ll2lL 12211 . It was predicted ll23ll that applying 



a small uniform bias current of positive or negative polarity, 
one can switch between the | J, and states, provided the dis- 
tance a between semifluxons is of the order of few (2 . . . 5) Aj. 
In this report, we demonstrate such a manipulation controlled 
by a uniform dc bias current and readout of two-semifluxon 
molecule states using on-chip SQUIDs. 

Experiment. Fig. [T] shows a sketch and a microscope im- 
age of one of the samples. To obtain a O-tt-O JJ we create 
two artificial phase discontinuities by using two pairs of cur- 
rent injectors attached to the top superconducting electrode of 
the JJ at a distance a from each other, see Fig. [TJb). The in- 
jector bias circuitry allows to feed both injector pairs by the 
current /jnj from a single current source, so that the phase 
discontinuities created by the injectors have opposite values, 
i.e. (k, — k), where k cx /jnj. Thus, one obtains a 0-k-Q JJ. 
As soon as K > 0, the Josephson junction reacts to the phase 
discontinuities by creating fractional vortices with topological 
charges that compensate these discontinuities, i.e. {—k,k). 
For < K < 2tt, one may have two possible stable antiferro- 
magnetically ordered fractional vortex molecules: a (— k, +k) 
molecule and a (2tt ~ k,k — 2tt) molecule. The first one has 
lower energy for k < tt, the second one for k > tt. For k ~ tt, 
these to states become degenerate ground states {—tt, tt) and 
(tt, — tt) of the system and are denoted as f j and J,|, respec- 
tively. The magnetic field oc d/i(a;) / dx in both states is shown 
in Fig. [Ha). The magnetic flux localized around each discon- 
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FIG. 1: (a) Sketch of the experimental circuit — a O-tt-O LJJ, posess- 
ing two degenerate ground states t| (red) and jt (black). The red 
and black curves were calculated numerically and represent the lo- 
cal magnetic field diJh{x)/dx. The magnetic field located at each 
phase discontinuity is coupled to and measured by a dc SQUID, (b) 
Optical microscope image of one of the investigated samples. Two 
current injector pairs are attached symmetrically to the top supercon- 
ducting electrode at a distance a (different from sample to sample). 
The widths Aw of each injector and the distance Aa- between them 
is Au) « Aa: ~ 2 ^ra. Two dc SQUIDs are placed in front of 
the phase discontinuities to read out the value of the flux localized 
at them. Each SQUID is coupled to the long JJ by using a super- 
conducting flux transformer (pickup loop in series with a single turn 
input coil). 



tinuity is coupled to one of the two washer-type dc SQUIDs 
by means of the pickup loop, see Fig.[TJb). The SQUIDs were 
designed to have a McCumber-Stewart parameter (3c ~ 1.1 
and normalized inductance parameter « 1.5. The work- 
ing point of each SQUID was chosen slightly above its critical 
current. 

The samples were fabricated at Hypres using a 30 A/cm^ 
Nb-AlOaj-Nb process. They all had L = 480 /im, but different 
values of a. Futher, the results are presented for the sample, 
which had a length L = 480 yum, and the distance between 
injectors a = 160 fim. The critical current measured T = 
4.2 K was 4 = 1.56 mA. 

First, to calibrate injectors, i.e., to find the value of /jnj 
for which k = tt, we have measured the /c(/inj) depen- 
dence, see Fig.|2a). The dependences Ic±{k) are 27r-periodic 
in K, so that from Fig. one can immediately say that 
^inj ~ 7.6 mA corresponds to k = 27t. The /c(/inj) curves 
were also calculated numerically for comparison (not shown). 
The most crucial parameter which is required to compare 
simulations with the experiment is the normalized JJ length 
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FIG. 2: (a) /c±(J'inj) dependences measured experimentally. The 
injector current of Jinj ~ 3.8 mA corresponds to the point k = tv. 
(b) The dependence of the SQUID voltage Vsq on Jinj measured 
experimentally at / = 0. 



I = L/Xj. We have estimated that I « 10.9 for Aj cal- 
culated using conventional expressions, and I w 7.5 for Aj 
calculated using idle region corrections 12411 . For our JJ ge- 
ometry the idle region corrections formula should give a bit 
overestimated values of I. Thus, the true value of I should 
be somewhere in between, but closer to 7.5. We have found 
that the simulated Ic (/i„j ) fits with the experimental one in the 
most reasonable way for L ~ 8.4Aj (a = 2.8A,/), which will 
be used below. In particular, for such a JJ length, the value 
of normalized critical cuiTent at k = tt of 0.58 is in excellent 
agreement with experiment, where 0.9/1.56 « 0.58. 

Second, without biasing the JJ (/ = 0), we have measured 
the left SQUID voltage T^L(/inj), see Fig.EJb). The VL(/i„j) 



and, therefore, $L(^inj), jump at certan values of /inj, corre- 
sponding to abrupt rearrangements of the flux in the molecule. 
By sweeping the /jnj in two directions, one can observe re- 
gions of bistability around k ~ ±7r. In these regions, shown 
by gold ribbons, both (— k, +k) and {—k + 2tt, k — 2tt) vortex 
molecules are stable, their energies and fluxes are comparable. 
The edges of bistable region seem to correspond to the edges 
of different branches of Ic± (/inj ) dependence. As we show 
later, it is just a coincidence. 

Next, we fix /jnj = 3.8 mA (k = it). A uniform bias cur- 
rent / exerts a Lorenz force acting on the fractional vortices 
and can induce transitions between the states || or |t if the 
force pushes vortices towards each other l23ll . Such a state 
transition takes place on the time scale of the inverse plasma 
frequency f^^ (in our case fp ^ 40 GHz) and will finally 
result in a reversal of the magnetic flux of both vortices in a 
molecule. While ramping the bias current / of the JJ back 
and forth, we measure the voltage V{I) across the J J as well 
as the voltages across each dc SQUID simultaneously, see 
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FIG. 3: (a) Voltage V{I) across the JJ at /inj = 3.8 mA (k = tt) 
and (b) corresponding dc SQUID voltages Vl{I) and Vr{I), mea- 
sured simultaneously. First, the bias current was swept with over- 
critical amplitude [light red and light blue reference curves in (b), 
black curve in (a)] and then with undercritical amplitude [red and 
blue curves in (b), not shown in (a) as V{I) = in this case]. The 
small interval of bias currents around J = for which the JJ was in 
the Meissner state is cut off from the reference curves in (b). Graphs 
(c) and (d) show recovered <I>'jJ'^(/) (red) and (blue) de- 

pendences and their comparison with numerically simulated <1>|J™(J) 
(gold) and $™(/) (green) curves. 



Fig.[3la,b). First, the bias current / was swept back and forth 
between the values exceeding the critical currents ±/c(/inj), 
so that most of the time the JJ was in the resistive (McCumber) 
state, see the black curve in Fig. [3 a) and light red and light 
blue curves in Fig. [3|b), further called "reference curves". 
Second, the bias current / was swept back and forth between 
the values not exceeding the positive and negative critical cur- 
rents Ic± (/j^ij ) of the JJ so that the JJ was in the zero voltage 
state during the whole sweep, see the red and the blue curves 
in Fig.[3|b). One can clearly see that the voltages across both 
S QUID s jump simultaneously at / w ±0.36 mA, indicating 
an abrupt rearrangement of the semifluxon molecule. 



In the resistive state, see the reference curves in Fig. Ob), 
the net magnetic flux coupled from the molecule to the 
SQUID averages to zero as a result of fast Josephson oscil- 
lations. Nevertheless, the reference curves in Fig. [3jb) show 
an almost linear Vl,r{I) dependence, because of the para- 
sitic coupling between the bias current / and the SQUIDs, 

'^LM Of I- 

The main result of the paper is presented in Fig. [3|c,d). It 
shows the magnetic fluxes sensed by the SQUIDs 

as a function of the (always undercritical) bias cuiTent /. 
The fluxes were calculated by subtracting the parasitic flux 
induced by /. To provide better fit, the flux axis was 
rescaled, which corresponds to the different coupling factors 
of SQUIDs. In these diagrams, two stable flux states with 
opposite fluxes and $'^^(7) are visible, correspond- 

ing to the states f | and it of the semifluxon vortex molecule. 
Transitions between these states were induced at the values 
of the applied bias current of I'^ w 0.36 mA. Note that 
the Josephson junction itself remains in the zero-voltage state 
during this measurement. The simulated depen- 
dences are also shown in Figs.[3];,d. The simulation was per- 
formed in normalized units for experimental parameters (L 
and a) using StkJJ.IZSI) The current axis was then rescaled 
taking into account the measured Ic of the sample. The re- 
arrangement current of /^g w 0.38 mA (in normalized units 
/^™//c = 0.25) obtained from simulation is in accord with 
I^^, see Figs.[3]:,d. 

In our experiment with electronically tunable discontinuity 
K oc /inj we have a unique possibility to study rearrangements 
not only in semifluxon AFM molecules (tt, — tt) ^ (— ti", tt), 
but also in molecules made of arbitrary fractional vortices 
(— K, +k) ^ (27r — K, +K — 27r), for a range of k in the 
vicinity of k = tt. To do this, a set of measurements similar 
to the one outlined above was performed for different values 
of K (X /inj . The values of bias current / at which the SQUID 
voltages Vl.b{I) jump were extracted from Vl,b,{I) curves 
using an automatic routine and plotted in Fig.|4|a) as a func- 
tion of /inj . The automatic routine was identifying the voltage 
jumps when a voltage change betwen neighboring points was 
exceeding a criterion of 0.4 /iV, chosen to be just above the 
voltage noise level. 

First, in Fig. Hla) one can see that both SQUIDs switch 
mostly simultaneously. In the regions where only the points 
corresponding to one SQUID are visible the other SQUID is 
insensitive to flux jumps and produces voltage jumps below 
the noise level. Second, one can observe several branches 
of the switching current vs. /inj. One branch corresponds 
to the switching at Ic and coinsides with the /^(/inj) curve 
also shown in Fig.Ua). Another two branches, parallel to the 
/inj axis, correspond to switching from the McCumber state 
to the superconducting state at return current / = i//?. Fi- 
nally, the two most interesting branches marked as /re (^inj ) 
look like a continuation of the /c(/inj) curves towards Ic = 0. 
They correspond to a rearrangement between the {k,—k) and 
(27r— K, +K— 27r) states. The magenta colored region between 
the /RE(-finj) lines represent a domain where the (k, — k) and 
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FIG. 4: (a) The Ic{Iinj) dependence measured using a voltage cri- 
terion algorithm in the vicinity of k = tt (black symbols). The val- 
ues of the switching current extracted from experimentally measured 
Vl(7) and Vr{I) curves taken at different linj are shown by red 
and blue symbols, (b) The critical current (black symbols) and re- 
arrangement currents (stars) calculated numerically for L = 8.4Aj 
and a = 2.8Aj. The region of coexistence of both states (bistability) 
(k, — k) and (k — 2n, —k + 2n) around k = tt is shown by magenta 
color. 



{2tt — K, +K — 2tt) states are both stable. It is this region, 
which can be used to build various switches, memory devices, 
bits and qubits. Fig.Ub) shows numerically calculated /c(k) 
and /re(k) dependences that show a state diagram very sim- 
ilar to the experimental one. 

Conclusions. We demonstrated experimentally that (a) a 
molecule of two interacting fractional Josephson vortices with 
two degenerate states || and |t can be created; (b) they 
can be conveniently manipulated between the f j and the || 
states by applying a small positive or negative dc bias cur- 
rent; and (c) the state of the molecule can be read out using 
on-chip SQUIDs. These results open a route to the applica- 
tion of fractional flux quanta as convenient information stor- 
age devices (superconducting memory). Simultaneously this 
demonstrates a ready to use preparation and readout technique 
for a semifluxon molecule based qubit i2U 12211 . 

Experimentally, one can easily fabricate molecules made 
of 3, 4, etc. fractional vortices or even a ID chain (crystal) 
made of many vortices. The state of these systems can also be 
controlled by applying a small bias current as presented here. 
Furthermore, the oscillatory eigenfrequencies of the vortices 
in the molecules or crystals can be controlled by the same bias 
current which in this case should be varied within the stability 
region of each state. Thus, a ID fractional vortex crystal can 



be used as a tunable band gap material lfl9ll . which opens the 
possibility to realize tunable filters in the sub-THz frequency 
range. 
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